Coreless Vortices in Rotating Two-Component Quantum Droplets 
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The rotation of a quantum liquid induces vortices to carry angular momentum. When the system 
is composed of multiple components that are distinguishable from each other, vortex cores in one 
component may be filled by particles of the other component, and coreless vortices form. Based on 
evidence from computational methods, here we show that the formation of coreless vortices occurs 
very similarly for repulsively interacting bosons and fermions, largely independent of the form of 
the particle interactions. We further address the connection to the Halperin wave functions of 
non-polarized quantum Hall states. 

PACS numbers: 67.10.-j 03.75.Lm 05.30.-d 73.21. La 



Setting a trapped Bose-Einstein condensate (BEC) ro- 
tating, as for example by stirring the bosonic droplet with 
lasers [1], quantum- mechanical "twisters" may form in 
the cloud to carry angular momentum. Such vortices 
consist of a hole in the particle density, associated with 
a quantized phase change in the order parameter around 
the core. With increasing rotation the vortices enter the 
bosonic cloud from its periphery [2|,|3[ one by one, finally 
leading to the well-known vortex lattice when angular 
momentum is high. The first experiment on vortices 
in a BEC [4] was in fact done with a two-component 
gas, as suggested by Williams and Holland [5[. Atoms 
of 87 Rubidium were trapped in two distinguishable hy- 
perfine spin states, and vortices could be created in one 
component while the other component was at rest. At 
high rotation, the vortex lattices of the two components 
are interlaced and 'coreless' vortices form with one com- 
ponent filling the vortex cores of the other component 
as a consequence of their mutual repulsion [6]. This was 
also experimentally observed by Schweikhard et al. [7|. 

While usually attributed to superfluidity, quantized 
vortices as they occur in condensates of bosonic atoms 
may even appear with trapped rotating fermion sys- 
tems [H, [lO] • Vortex patterns in few-electron quantum 
dots at strong magnetic fields appear as hole- li ke q uasi- 
particles in much analogy to the bosonic case 
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Moreover, there often exists a simple transformation of 
the fermionic wave function into the bosonic one [l2[ and 
a close theoretical connection between rotating bosonic 
systems and quantum Hall states [l3[. Depending on 
the angular momentum vortices may either nucleate to 
a vortex cluster, just like in the bosonic case @, 0, flij . 
or appear as particle- vortex composites in the finite-size 
precursors of polarized [1, H, [Io|, EH and unpolar- 



ized [la, ll7[ fractional quantum Hall states. In experi- 



ments rotating two-component fermion droplets may be 
realized with trapped dipolar fermion gases [l8[ , or with 
electrons in quasi two-dimensional quantum dots [19[, 
where the two components are distinguished by the spin 



quantization axis of the electrons, and coreless vortices 
are expected to form in the regime of rapid rotation. 

In this Letter, we report numerical evidence for the 
remarkable similarity of coreless vortices in rotating two- 
component quantum droplets regardless of the particle 
statistics and the form of particle interactions. 

We focus on rapid rotation, and restrict our analy- 
sis to a two-dimensional harmonic trapping potential of 
oscillator frequency uo. For a droplet composed of two 
distinguishable species A and B with particle number 
N — Na + Nb , the Hamiltonian is 
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where r = (x,y), m is the mass (equal for both species), 
Q is the frequency of the external rotation, and L is 
the (z-component) of angular momentum. V is the 
particle-particle interaction potential, which is Coulom- 
bic, V = e 2 /(47re|r^ — rj|), in the case of electrons in 
a quantum dot. For bosonic condensates we have used 
both the Coulombic potential and the contact potential 
V = U5(yi — rj), where U = Air^a/m and a is the 
scattering length (here assumed to be the same for both 
components), which is more realistic for cold atom con- 
densates. In semiconductor quantum dots, the Zeeman 
splitting can be engineered to be zero, as for example in 
GaAs/Al^Gai-^ As-based systems [20I, [2l[. I n this paper 
we explicitly assume that the Zeeman splitting is neg- 
ligible which leads to symmetry between the two spin 
components. 

The eigenstates of the Hamiltonian Eq. (pQ) are calcu- 
lated by numerical diagonalization in the lowest Landau 
Level (LLL) (for details see Refs. [loj and [22[). For large 
systems, however, we applied the spin-density-functional 
method as described in Ref. [8[. 

Pair-correlation functions are commonly used to ana- 
lyze the structure of the many-body wave function but 
integrations average out the effect of phase singularities 
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FIG. 1: ( Color online): Identification of vortex structures in a 
restricted wave function. The figures show RWFs for spin- 1/2 
fermions. The triangles show the fixed particle positions and 
the probing particle is marked with light blue. Contours show 
the magnitude of the RWF (on a logarithmic scale), and the 
gray- scale shows the phase (darkest gray = 0, lightest gray = 
2tt). a) Vortex which is not attached to any particle, b) Pauli 
vortex (blue) which is mandated by the wave function anti- 
symmetry between interchange of indistinguishable fermions. 
c) Coreless vortex where a vortex core in one component is 
filled by a particle in the other component, d) A particle- 
vortex composite made up of a fermion (with a Pauli vortex) 
and two additional nodes which are bound to the particle. 



which are signatures of vortices. Therefore, we investi- 
gate the nodal structure of the wave function using re- 
stricted (conditional) wave functions (RWF) [8] which 
essentially depict the structure of the 27V-dimensional 
eigenstate of the Hamiltonian Eq. (pQ), \£(ri, . . . , r/v), in 
a two-dimensional subspace by fixing all but one probing 
particle: 
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Here, the denominator is for normalization and the fixed 
particle positions , . . . , are determined such that the 
norm of ^RWF is maximized and is the most prob- 
able position of the probing particle. The RWFs are 
two-dimensional complex functions where vortices can be 
identified directly as nodes with a phase change of a mul- 
tiple of 2tt in a path around them. Figure [T] provides a 
guide to the identification of vortex structures in RWFs. 

In the following, we exclusively discuss equal pop- 
ulations of the two species A and B. We first study 
coreless vortex structures in a 6-particle system with 
Na — Nb = 3, using the ED method in the LLL. 
For a single-component fermion system the state with 
the lowest angular momentum in the LLL is called the 
maximum-density-droplet (MDD) [23], and has angular 
momentum Lmdd = N(N — l)/2, which is the mini- 
mum angular momentum compatible with the Pauli prin- 



ciple. The corresponding wave function is ^mdd = 
^-i<j( z i ~ z j) e ~ (with z = x + iy), where the zeros 

in the polynomial are interpreted as Pauli vortices [§]. 

For single-component systems a simple transformation 
removes the attachment of Pauli vortices to each fermion 
and thus leads to a symmetric wave function [24[ which 
usually describes, to a significant degree of accuracy [l2| . 
bosonic states at Lbosoii = £ Fermion _ ^mdd- In two- 
component systems Pauli vortices are only attached to 
particles of the same spin and there exist states without 
bosonic counterparts at £ Fermion < ^mdd- Thus, we 
consider only the angular momentum regime ^Fermion > 
£mdd- 

Figure [2] shows the RWFs of two-component bosonic 
and fermionic systems for Coulomb interactions, at 
ground-state angular momenta as given in each panel. 
Removal of a single vortex from each fermion position (for 
both components) leaves a wave function whose nodal 
structure shows formation of an equal number of coreless 
vortices in both fermionic and bosonic systems. These 
can be identified in the RWFs as vortices in one compo- 
nent that coincide with probability maxima in the op- 
posite component, in order to minimize the interaction 
energy. States at intermediate angular momenta show a 
vortex entry from the drop let p eriphery, similar to that in 
single-component systems [l5ll25j|. Just as in the bosonic 
case, the two components of a fermion system first sep- 
arate (Fig. [2)3) [22[ and, as the angular momentum in- 
creases, form interlaced vortex patterns (Fig. [2]i). These 
features in the RWFs are well reflected also in the pair 
correlations of bosonic systems (Fig. [3]). Furthermore, 
we find that for bosons interacting with a contact inter- 
action the overlaps with the corresponding states with 
Coulomb interaction are higher than 0.98 which shows 
that the coreless vortex formation is largely independent 
of the form of the interaction. 

There is a trend for the localization of particles to in- 
crease with angular momentum (and thus vorticity), as 
manifested by a narrowing of the RWF peak around the 
most probable particle position (Fig. [2|). The effect is 
more clearly seen in fermion systems where the number 
of vortices is higher. 

For larger system sizes, the diagonalization of the 
many-body Hamiltonian Eq. (pp) becomes numerically 
impossible. In bosonic systems one often turns to 
mean-field methods such as the Gross-Pitaevskii (GP) 
equation. In the case of fermions, we can make use 
of density- functional theory (DFT), as previously de- 
scribed for one-component droplets in Ref. [8]. A single- 
determinant wave function constructed from the Kohn- 
Sham orbitals [26[ shows formation of coreless vortices 
in much analogy to the ED results in Fig. [2Jd. Figure 
[H shows the separation of the two components with a 
sharp phase boundary, analogously to the bosonic case. 
For states with high angular momentum, vorticity be- 
comes large and the DFT solutions show interwoven vor- 
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FIG. 2: (Color online): The restricted wave functions calcu- 
lated with the exact diagonalization method for a two- com- 
ponent system with Na = Nb = 3 for bosons at angular mo- 
menta Lb (left) and fermionic systems at angular momenta 
Lf = Lb + ^mdd = Lb + 15 (right). The states shown are 
ground states in the lowest Landau level approximation. The 
bosonic and fermionic states are analogous in nodal structure 
when one node is removed from each fermion (these nodes are 
not marked in the right panel for clarity). The figure shows 
the attachment of b) one, c) two, and d) three vortices, re- 
spectively, to particles of the opposite species which indicates 
the formation of coreless vortices. The notation is explained 
in the caption of Fig. [T] 
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FIG. 3: (Color online): Difference of pair-correlation func- 
tions for two species of bosons, for Na = Nb = 3. a) Com- 
ponents are separated due to formation of the first coreless 
vortex, b) Another domain forms due to formation of the 
second coreless vortex, c) Formation of the third coreless vor- 
tex is associated here with "antiferromagnetic" ordering as 
seen in Fig.EJl. The reduced wave functions of corresponding 
states with similar symmetries are shown in Fig. [2j>d. 



N=20 



Spin density 





FIG. 4: (Color online): The DFT spin-density (left) in a 20- 
electron droplet with equal number of spins (iVj = N± = 10) 
shows separation of spin components into domains. The an- 
gular momentum L — 193 is slightly above the MDD angular 
momentum Lmdd = 190. The restricted wave function for 
the spin-down electrons (right) shows a nodal structure with 
nodes equal to the number of fixed particles and an additional 
coreless vortex entering the droplet (arrow). 



tex sheets as previously reported for bosonic systems us- 
ing the GP equations [27[ . 

In the regime where the number of vortices in a bosonic 
system exceeds the number of particles, the system shows 
a mixture of coreless vortices and particle-vortex com- 
posites. Particularly, the bosonic ground state at L = 21 
shows composites of particles with two vortices in the 
component of the probing particle, and coreless vortices 
in the other component, as seen in Fig. [Sfci- 

The structure of the many-body states analyzed here 
bear much resemblance with the trial wave functions of 
non-polarized quantum Hall states. Quantum Hall states 
are usually polarized in high magnetic fields due to the 
Zeeman-coupling but some quantum Hall plateaus, such 
as v = | and v = |, have been proposed to be non- 
polarized states 12 8l 1 291 . The Halperin wave function for 
such states reads 1 301 
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where q is an odd integer (due to fermion antisymmetry) 
and p is a positive integer and the Gaussians have been 
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FIG. 5: (Color online): a) In the regime of rapid rotation the 
bosonic ground state at L = 21 shows a mixture of both core- 
less vortices and particle-double- vortex composites b) L = 24 
fermion state with a nodal structure which corresponds to 
that of the Halperin-wave function with p = 2, q = 1, which 
attaches two vortices on top of electrons of the opposite spin 
and one (Pauli) vortex on top of electrons of the same spin. 
Note that in contrast to Fig. [2] we here mark all nodes. 

omitted. This wave function attaches different number of 
vortices on top of spin up and spin down electrons. A vor- 
tex attached to the opposite spin component can be inter- 
preted as a coreless vortex. We find that this nodal struc- 
ture has a finite-size analogue in two-component fermion 
droplets. In Ref. [22[ the L = 24 state at N = 6 was inter- 
preted as a finite-size counterpart of the v = | quantum 
Hall state with a Halperin wave function p = 2, q = 1. 
Figure [5)3 shows the RWF of this state where one (Pauli) 
vortex is attached to each particle of the same spin and 
two (coreless) vortices are attached to particles of the op- 
posite spin, in good agreement with the Halperin model. 
However, the overlap of this state with the Halperin wave 
function has been found to be low for large particle num- 
bers due to mixing of spin states in the Halperin wave 
function (Ref. [U). 

We finally note that even though this study focuses on 
systems with equal component sizes, coreless vortices also 
form with unequal component populations. In bosonic 
systems these have been studied extensively m. As a 
fermionic example we consider results in Ref. [3l| where 
the MDD breakup was analyzed in quantum dots and 
some beyond-MDD states were found to have partial spin 
polarization. On the basis of the results presented above, 
these states can be re-interpreted as vortex states in the 
majority spin component whose vortex cores are filled 
with an electron of the minority spin component. Spin- 
selective electron transport experiments could be able to 
indirectly probe for these types of coreless vortex states 
in quantum dots [32[. 
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